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A mixed supersymmetric-algebraic approach to construction of the minimum uncertainty coherent 
states of anharmonic osciUators is presented. It permits generating not only the weU-known coherent 
states of the harmonic and Morse osciUators but also the so far unknown coherent states of the 
Wei Hua, Kratzer-Fues and generalized Morse and Kratzer-Fues oscillators. The method can be 
applied also to generate superpotentials indispensable for deriving the Schrodinger equation in the 
supersymmetric form amenable to direct solution in the SUSYQM scheme. 
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c : 

^ . The coherent states introduced by Schrodinger in 1926 [1], make a very useful tool for investigation of various 
' ' problems in quantum optics [2] , in particular the interactions of matter with coherent radiation [3] . Such resonant 
interactions produce the coherent effects, e.g. self-induced transparency, soliton formation, excitation of a coherent 
' superposition of rotational states, periodic alternations of the refractive index in both nuclear and molecular systems 
[3]. Studies in the latter case require construction of the coherent states for anharmonic molecular potentials. Such 

■ states are defined in the similar manner as the ordinary coherent states of the harmonic oscillator [4]: (i) they are 
eigenstates of the annihilation operator, (ii) they minimize the generalized position-momentum uncertainty relation 

-ilj • and (iii) they arise from the operation of a unitary displacement operator to the ground state of the oscillator. It 
C ' should be pointed out that the definition (iii) relies on the form of the displacement operator, which is specific to 
2 I harmonic oscillator [5], hence, in this case mainly approximate coherent states can be derived using, for example, 
qh' Nieto and Simons [6] or Kais and Levine [7] procedures. 

I— I , Coherent states of anharmonic oscillators have been constructed using several alternative approaches. In the method 
proposed by Nieto and Simmons [6] , the position and momentum operators are chosen in such a way that the resultant 
' ^ Hamiltonian resembles that for an harmonic oscillator. The coherent states are then determined on condition that 
^ I . they minimize the generalized uncertainty relation in the new variables. Perelomov [8] has derived the coherent 
ly-j [ states using the irreducible representations of a Lie group. This method has been successfully applied to generate the 
QQ : coherent states of the Morse oscillator [7] . The generalized coherent states can also be generated using an algebraic 
CO ' method [5] and supersymmetric quantum mechanics (SUSYQM) [9]. Applying the above formalisms, the coherent 
\^ , states for Morse [5, 6, 7], Poschl- Teller [11], hydrogen atom [12], Eckart and Rosen-Morse [13], double- well and linear 

• (gravitational) potentials [14] have been constructed. 
[ — . , In the present study, we introduce the mixed supersymmetric-algebraic method, which permits generating not only 

■ the coherent states of the harmonic and Morse oscillators but also the so far unknown coherent states of the Wei Hua 
[15], Kratzer-Fues [16] and generalized Morse and Kratzer-Fues oscillators. The method starts from the vibrational 
dimensionless Schrodinger equation 



X 



If + Viq) - Eo 



\v>=AEyo\v>, ^ = 



in which AE^q = Ey — Eq whereas q — u^r denotes a dimensionless spatial variable r, with a scalling factor 
depending on the explicit form of the potential energy term V{q). 

The crucial for the approach proposed is the assumption that the last two terms in the operator part of Eq. (1) 
can be specified in the form of the Riccati equation 



V[q) 'E^ = \ 



dq 



(2) 



familiar in SUSYQM [17]. Here, x{q) is the anharmonic coordinate, which satisfies the commutation relation [x{q),p\ = 
idx{q) / dq. Its form depends on the oscillator type, hence the explicit expression for x{q) will be determined for a given 
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form of the potential energy function. In SUSYQM x{q) (with accuracy to sign) is interpreted as a superpotential 
[17], which permits construction of the supersymmetric Schrodinger equation straightforward to analytical solutions. 
Substituting Eq. (2) into (1) one gets the latter in the factorized form 



A^A\v >= AEyo\v >, 



in which 



V2 



\/2 



A, A'' 



dx{q) 
dq 



(3) 



(4) 



In order to construct the coherent state for the potential V{q), we need a ground state solution |0 > of Eq. (3), which 
is an eigenstate of the operator A. If A annihilates the ground state A\0 >= 0, then the coherent states \a > are the 
eigenstates of the annihilation operator A and the following relations are fulfilled 

A\a >= a\a >, < a\a* =< a\A^, \a >= |0 > exp(\/2ag). (5) 
One may prove that the states \a > minimize the generalized position- momentum uncertainty relation [5] 



{Ax{q)fiApf > -< a\g{q)\a >^ g{x) = -i[x{q),p] 



dx{q) 
dq 



To prove this thesis let's calculate 



< a\x{q)\a >= < a\A + 4^|q! >= ;^(" + a*), 



(6) 



(7) 



< a\p\a >= < ~ >= ~ 



(8) 



2 < a\x{qf\oi >=< a\AA + lA^ A + A^ A^ 



dx(q) 
dq 



a >= 



(a + a*)^- < a 



dx{q) 
dq 



a > 



(9) 



- 2 < a\f\a >=< alAA - 2A^A + i^it + '^^\a > 

dq 

in which Eq. (4) is employed in the form AA^ = A — dx{q)/dq. 
Having calculated (7) - (10), we can pass to evaluate 



(a - a*y+ < a 



dx{q) 
dq 



a > 



{Ax{q)f =< a\x{q)f\a > - < a\x{q)\a >^= ~ < a|^^^|a >, 



(10) 



(11) 



(Ap)^ =< a\p^\a > - < a\p\a >^= -- <a 



dx{q) 
dq 



a > 



providing Ax{q) = Ap and 



{Ax{q)f{Apf = \<a 



I dx{q) 
' dq 



a 



(12) 



(13) 



The calculations performed prove that the states \a > minimize the generalized position- momentum uncertainty 
relation for the anharmonic coordinate x{q). They are also the eigenstates of the operator A, which annihilates the 
ground state A\0 >= 0, hence they satisfy the two fundamental requirements established for the coherent states of 
an anharmonic oscillator. 

In order to demonstrate how the method works, let's calculate first the well-known coherent states of the harmonic 
oscillator. To this purpose let's assume that 



dx(q) 
dq 



-l^x{q) = -q for a;(0) = 0. 



(14) 
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Then Eq.(l) including (2) turns out to be the well-known Schrodingcr equation for the ground state |0 >= cxp(— g^/2) 
of an harmonic oscillator, whereas the operators A and ^4^ take the well-known form of annihilation and creation 



operators, which satisfy the commutation rule 



oscillator can be specified by the general formula (5) 



-dx{q)/dq = 1. Hence, the coherent states of an harmonic 



\a >= |0 > exp(-\/2ag) = exp(— ig^) exp(-\/2Q;g). 



(15) 



The results obtained indicate that crucial for the method proposed is the explicit form of the term dx{q)/dq, which 
in the general case can be given as a negative x-dependent function 



dx{q) 
dq 



(16) 



Hence, employing different analytical functions f{x) one may generate a variety of potentials and associated super- 
potentials satisfying relation (2). The indispensable for this purpose x{q) can be calculated from (16) by integration, 
provided that we known the explicit form of f{x). Assuming that \x{q)\ < 1 one may expand f{x) in a power series 



f{x) = ci{x + co/ci) -I- C2{X + co/cif + . 



(17) 



and then successively apply the first-, second- and higher-order terms in determination of coherent states of anharmonic 
potentials. For example, if we take into account the linear term f{x) = ci (a; -|- cq /ci) and assuming x{0) = (1 — cq) / ci 
one may calculate from (16) 



x{q) = — [exp(-cig) - Co] . 

^1 



(18) 



Introducing (18) into (2) and redefining the constants cq = s — Xe, c\ = \/2xl from Eqs. (1) and (2) one gets the 
Schrodinger equation 



1 ( d^ 1 r , ,t2 a^el ,„ 

2|-d?+2^[^-^^P(-^^)] -a+-||0>=0 



and the ground state solution 



|0 >= exp[— - — exp{—^/2xeq)] exp 

ZXfi 



/2x. 



:(S - Xe) 



(19) 



(20) 



of the generalized Morse oscillator [10] V{r) = De[s — exp(— ar)] . For s = 1 the above specified formulae are reduced 
to the well-known equations derived by Cooper [5], in which Xe = JiUe/4:De denotes the anharmonicity constant, 
We = ay^2De/m is the vibrational frequency defined by the reduced mass m of the system, the dissociation energy 
Df. and the range parameter a appearing in the Morse potential. 

Eqs. (18) and (4) produce the generalized Morse annihilation and creation operators 



V2 



d 
dq 



and associated coherent states 

\a >= exp 



'QiXf' 



-^exp(-V2xIg) 



' V2 



d 
dq 



2iCe 



exp 



fix. 



(s - x^) 



exp 



(\/2Q;g). 



(21) 



(22) 



They are eigenstates of the annihilation operator ^4, which minimize the uncertainty relation (6) for 

exp(-A/2i7g)- 

Taking into account the parabolic expansion 



f{x) = ci{x + Co/Ci) + C2{X + co/cif 



(23) 
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and the identical assumption as previously x(0) = (1 — co)/ci, from Eq.(16) one obtains 

(cci/c2)exp[-ci(g- go)] cq 



x{q) 



1 - cexp[-ci(g- qo) 



Cl 



(24) 



in which c = C/{B/W-C), go - ln(B/VK-C)/ci, = (2co + c?)/[2ci(l-C2)], B = Ci/(c? + C2) and C = C2/(cf +C2). 
Hence, the Schrodinger equation (3) and its ground state solution take the forms 



2D 



l-exp[-ci((?-(?o)] 



2Eq ) |0 >= 0, 



1 - cexp[-ci(q - go) 

-L. ^ 

>= {1 - cexp[-ci(g - go)]}'^^ {ccxp[-ci(g - go)]}"i , 



(25) 



(26) 



in which 2D = (1 — C2)W^, 2Eq = — c^/c\. Eq. (25) is a well-known quantal equation for the ground state of the 
Wei Hua oscillator [15], whose coherent states have not been derived as yet. Applying the above specified procedure 
one gets the correspondig annihilation and creation operators 



\/2 



it 1 
\/2 



d_ _ (CC1/C2) exp[-ci(g - go)] _^ Cq 
dq 1 - cexp[-ci(g - go)] ci 



_ (CC1/C2) exp[-ci(g - go)] ^ Cq' 
dq 1 - cexp[-ci(g - go)] ci 



as well as the coherent states of the Wei Hua oscillator 

\a {1 - ccxp[-ci(g- go)]}~ {ccxp[-ci(g - go)]}"i cxp(\/2Q:g). 



(27) 
(28) 

(29) 



They are cigcnstatcs of the annihilation operator A\a >= a\a > and minimize the uncertainty relation (6) for 

A, = {ccl/c2) exp[-ci(g - go)]/{l - cexp[-ci(g - go)]}^- 

The method proposed also permits a derivation of the coherent states of the Kratzer-Fues potential [16] specified 
in the form 



V{r) = D, 



> - re' 


1 








z 


r 


A + z. 



(30) 



in which denotes the equilibrium configuration of the molecule bond V{re) = 0, whereas z = (r — r^jr,, is the 
Dunham variable. To this aim we employ the generating function j{x) = [ci{x + co/ci)]^, which introduced into Eq. 
(16) yields 



a;(g) 



Co 



x(0) = (l-co)/ci. 



ci(cig-l-l) Cl' 

Then the Schrodinger equation (1) and its ground state solution for x{q) specified above take the form 



cig - s 



1 + cig 



2Eo \ |0 >= 0, |0 >= (1 + cig)^ exp 



(l-c2)(cig + l) 



c?(s + l) 



(31) 



(32) 



in which 2D = Cg/fcf (1 - cf)], s = {I-cq- cD/cq, 2Eo = Cq/{1 - cf). It is easy to verify that for s = or cq = 1 - cf 
and ciq = z, Eq. (32) turns out to be the Schrodinger equation for the Kratzer-Fues oscillator [16] 



2 I rfg2 



2D 



Cig 



1 2 



1 + cig 



2Eo}\0>=0, |0>= (l + cig)=i exp[-(l - c2)(cig + l)/c?]. 



(33) 



in which 2D = (1 — cD/ci, 2Eq = 1 — cl. Hence the Kratzer-Fues annihilation and creation operators can be given 
in the form 



V2 



1 



dq ci(cig-l-l) Ci 



it- — 
' V2 



1 



dq ci(cig-l-l) Ci 



(34) 
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and the associated coherent states 



\a >= (1 + cig)°? exp[-(l - cf)(cig + l)/cf\ exp{\/2aq) 



(35) 



minimize the uncertainty relation (6) for 



^ =(cig + l)-2. 

The potential V{z) = D,,[{z — s)/(l + z)]"^, which appears in (32) is worth considering as it represents a generalized 
version of the Kratzer-Fues formula V{r) = D(,[l — re{s + l)/r]^. From the relation — 1 < 1 — re(s + l)/r < 1 one 
may calculate the convergence radius R G [re{s + l)/2, oo] for the new potential. The former increases for s G (—1, 0) 
in comparison to the original Kratzcr-Fues potential (s = 0) yielding R S (re/2,oo). In such circumstances, the 
expansion of the potential energy function 



V{r) = CO 



.(s + 1) 



N 
n=l 



.(s + l) 



(36) 



into a series of the generalized Kratzer-Fues variable, will provide much accurate reproduction of the real potential 
curves then that obtained by the Simons-Parr-Finlan expansion (s = ) [18], which diverges in the imited-atom limit 
r — » 0. The set of parameters (re, s, cq, ci, ...) can be evaluated from the molecular IR and MW spectra by making 
use of the fitting procedure. It should be pointed out also that the new potential permits analytical solution of the 
Schrodinger equation and can be used to generate the coherent states for the generalized Kratzer-Fues oscillator 



\a >= (1 + cig) °? exp[-(l - cl){ciq + l)/c?(s + 1)] exp{V2aq). 



(37) 



The method proposed is general and permits construction of the coherent states, associated potentials and superpo- 
tentials as well as deriving the supersymmetric Schrodinger equation amenable to direct solution in the SUSYQM 
scheme [19]. In the standard approach the superpotentials are solutions of the Riccati equation (2) obtained for 
the specific form of the potential function V{q) [19]. Here a new procedure has been introduced, which permits 
simultaneous derivation of the potentials and associated superpotentials assuming that dx{q)/dq = —f{x). This term 
can be expanded in a power series of x{q) and then used to generate the coherent states for different orders of the 
expansion (17). For the terms up to the second order, the method produces the minimum uncertainty coherent states 
for harmonic, Morse, Wei Hua, Kratzer-Fues and generalized Morse and Kratzer-Fues potentials. They are the most 
important potential energy functions employed in molecular quantum theory, theoretical spectroscopy and quantum 
optics. 



[1] E. Schrodinger, Naturwiss. 14, 664 (1926). 

[2] R. J. Glauber, Phys. Rev. Lett. 10, 84 (1963); Phys. Rev. 130, 2529 (1963); 131, 2766 (1963). 

[3] J. P. Heritage, T. K. Gustafson and C. H. Lin, Phys. Rev. Lett. 34, 1299 (1975); R. B. Walker and R. K. Preston, J. 
Chem. Phys. 67, 2017 (1977); N. Mankoc-Borstnik, L. Fonda and B. Borstnik, Phys. Rev. A 35, 4132 (1987); L. Fonda, 
N. Mankoc-Borstnik and M. Rosina, Phys. Rep. 158, 159 (1988); N. M. Avram, Gh. E. Draganescu and C. N. Avram, J. 
Opt. B: Quantum Semiclass. Opt. 2, 214 (2000). 

[4] W-M. Zhang, D. H. Feng and R. Gilmore, Rev. Mod. Phys. 62, 867 (1990). 

[5] I. L. Cooper, J. Phys. A; Math. Gen. 25, 1671 (1992). 

[6] M. M. Nieto and L. M. Simmons, Phys. Rev. Lett. 41, 207 (1978); Phys. Rev. A 19, 438 (1979). 

[7] S. Kais and R. D. Levine, Phys. Rev. A 41, 2301 (1990), G. C. Gerry, Phys. Rev. A 31, 2721 (1985), G. C. Gerry, Phys. 
Rev. A 33, 2207 (1986). 

[8] A. M. Perelomov, Generalized coherent states and their applications (Berlin, Springer, 1986). 

[9] T. Pukui and A. Aizawa, Phys. Lett. A 180, 308 (1993); M. G. Benedict and B. Molnar, Phys. Rev. A 60, R1737 (1999). 

[10] P. M. Morse, Phys. Rev. 34, 57 (1929). 

[11] M. M. Nicto and L. M. Simmons, Phys. Rev. D 20, 1332 (1979); M. G. A. Crawford and E. R. Vrscay, Phys. Rev. A 57, 

106 (1998) ; A. H. Kinani and M. Daoud, Piiys. Lett. A 283, 291 (2001). 
[12] C. C. Gerry and J. Kiefer, Phys. Rev. A 37, 665 (1988); J. R. Klauder, J. Phys. A: Math. Gen. 29, L293 (1996); P. 

Majumdar and H. S. Sharatchandra, Phys. Rev. A 56, R3322 (1997); S. A. Pol'shin, J. Phys. A: Math. Gen. 33, L357 

(2000). 

[13] A. Chenaghou and H. Fakhri, Mod. Phys. Lett. A 17, 1701 (2002). 
[14] M. M. Nieto, Mod. Pliys. Lett. A 16, 2305 (2001). 
[15] W. Hua, Phys. Rev. A 42, 2524 (1990). 

[16] A. Kratzer, Z. Phys. 3, 289 (1920); E. Fues, Ann. Phys. Paris, 80, 281 (1926). 
[17] F. Cooper, A. Khare and U. Sukhatme, Phys. Rep. 251, 267 (1955). 
[18] G. Simons, R. G. Parr and J. M. Finlan, J. Chem. Phys. 59, 3229 (1973). 
[19] F. Schwabl, Quantum Mechanics (Berlin, Springer, 1992). 



